Dependence of hadronic properties on Quark Masses and Constraints on their 

Cosmological Variation 
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We follow our previous paper on possible cosmological variation of weak scale (quark masses) 
and strong scale, inspired by data on cosmological variation of the electromagnetic fine structure 
constant from distant quasar (QSO) absorption spectra. In this work we identify the strange quark 
mass m s as the most important quantity, and the sigma meson mass as the ingredient of the nuclear 
forces most sensitive to it. As a result, we claim significantly stronger limits on ratio of weak/strong 
scale (W = m a / Nqcd) variation following from our previous discussion of primordial Big-Bang 
Nucleosynthesis (\6W/W\ < 0.006) and Oklo natural nuclear reactor (\6W/W\ < 1.2 ■ 10~ 10 ; there 
is also a non-zero solution SW/W = (-0.56 ± 0.05) • 10" 9 ) . 



I. INTRODUCTION 

The understanding of how exactly the fundamental 
parameters of the Standard Model enter any observ- 
able is certainly one of the most important aims of 
hadronic/nuclear physics. Two old deep questions drives 
its discussion: (i) Can there be "alternative Universes" 
with different set of parameters, and what are the bound- 
aries of the World we know in the parameter space?; 
(ii) How to observe cosmological variations of weak and 
strong scales? 

Discussion of both issues has been significantly revived 
recently. We will not go into questions (i) (see e.g. [1]) 
and only mention the latter (ii). The issue of cosmo- 
logical time variation of major constants of physics has 
been recently revived by astronomical data which seem 
to suggest a variation of electromagnetic a at the 10~ 5 
level for the time scale 10 bn years, see [2]. The statis- 
tical significance of the effect at the moment obviously 
exclude any random fluctuations, so the effect definitely 
exists. Whether it may or may not have a conventional 
explanation is not yet clear: more experimental work is 
clearly needed to reach any conclusions. 

Nevertheless, it is quite timely to have another look at 
existing limits on time variation of all the fundamental 
constants. In particular, since the electromagnetic and 
weak forces are mixed together in the Standard Model, 
one may expect a similar modification of the weak cou- 
plings, the weak scale in general and quark masses in 
particular. In fact, one can measure only variation of 
dimensionless parameters. Therefore, we obtain limits 
on variation of m s j Aqcd where m s is the strange quark 
mass and Aqc d is the QCD scale defined as a position of 
Landau pole in logarithm for running coupling constant. 
It is convenient to put Aqcd = const. 

A generic further argument goes as follows. The 
masses of 3 heavy quarks - c,b,t - are too large to be 
important in hadronic and nuclear physics. The masses 
of two light ones - u,d - are important, in particular 
via the pion contribution to nuclear forces we studied in 



our previous paper [3] and a subsequent one [4]. Much 
more extensive discussion of this issue in a contexts of 
chiral perturbation theory can be found in literature, see 
e.g. [5] and references therein. The conclusion reached 
in those studies is that the m u ,md are too small to be 
really important. 

Thus, we focus on the dependence on the strange 
quark mass m s , the only one which has the right 
magnitude to generate a maximal sensitivity of the 
hadronic/nuclear physics to the weak scale. Indeed, as 
follows from QCD phenomenology, its variation from 
to experimental value influence vacuum parameters such 
as the quark condensates qq at the factor-2 level. It 
also affects the masses of even non — strange hadrons 
such as the nucleon at about 20 percent level, etc. 
The fundamental explanation of why such unexpected 
"strangeness" is in fact present is related with the impor- 
tant role of the instanton-induced effects in QCD. As it 
follows from the chiral anomaly relation and is explained 
in details e.g. in a review [6], the multi-fermion 't Hooft 
interaction necessarily involves all 3 flavors at all steps, 
even in the interactions of light u,d quarks. 

The main new element of the present work is iden- 
tification of the specific hadron, the cr, now known as 
/o(600) meson in the Review of Particle Properties [10], 
as the important ingredient of the nuclear physics most 
sensitive to m s . This resonance is known for a very long 
time as a structure in irir scattering, and is widely used 
in nuclear physics, see e.g. a review on applications of 
Walecka model [7]. And nevertheless, its acceptance as a 
real hadronic state was difficult and had a very compli- 
cated history. Non-relativistic quark models of hadrons 
intensify it is as a 1=1 or p-state, and tend to predict 
its mass to be around 1.4 GeV or so, in striking contra- 
diction to experiment which finds a mass of only about 
m a ~ 500 MeV. Also large width of a makes this state 
to be easily deformed by all kind of effects, see e.g. re- 
cent work by one of us predicting drastic modification on 
the sigma shape in pp and heavy ion collisions [9]. On 
the other hand, recent data have elucidated sigma pro- 
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duction in a set of much simpler situations, from heavy 
quark hadrons, T transitions and D decays. Those con- 
sistently point toward the smaller width r CT ~ 250 MeV . 
This reduced the controversy and put this particle back 
into the Review of Particle Properties [10]. 

In the next section we will argue that the sigma mass 
has stronger sensitivity to m s , than that of ordinary non- 
strange hadrons like N or u). This happens because of a 
valence strange part plus the repulsion from the nearby 
KK,r\r\ continuum thresholds. We then estimate the 
derivative of the deuteron binding and neutron resonance 
energies to m s . 

II. SENSITIVITY OF HADRONIC MASSES TO 
VARIATION OF QUARK MASSES 

A. Why sigma? 



The very important lesson about nuclear forces this 
model had emphasized is that the nuclear potential is 
in fact a highly tuned small difference of two large terms. 

We will argue in the next section that there are reasons 
to think that the sensitivity of these two terms to the fun- 
damental weak scale is quite different. Sigma (scalars) 
involve all quark flavors while omega (vectors) do not, 
forbidden by Zweig rule. As a result, there is no fine 
tuning in the derivative over m s , which significantly en- 
hances the effect to be derived. 

The values of the two coupling constants used in the 
nuclear matter applications of this model [7] are 

g\ = 357Aml/m 2 N w 100 

gl = 273.8ml/m 2 N w 190 (3) 

for m a = 500 MeV. Note that the effective scalar cou- 
pling is close to naive value (1) mentioned above. 



The very first appearance of the sigma mesons was 
as a two-pion scalar-isoscalar resonance. It has been 
gradually learned that the corresponding channel for qq 
interaction is the "maximally attractive channel" , with 
attraction so strong that it breaks spontaneously chiral 
symmetry and produce the non-zero quark condensate. 
The mechanism of that attraction is attributed mostly to 
instanton-induced 't Hooft interaction, for a review see 
[6]._ 

Sigma meson is an excitation on top of the scalar con- 
densate, a kind of a Higgs boson of strong interactions. If 
indeed one naively assumes that it underlines all hadronic 
masses, e.g. that of the nucleon, the corresponding cou- 
pling can be estimated as 



g s = M n /U « 10 



(1) 



. This large value in turn implies that the perturbation 
theory can only be used as a qualitative guide, at best. 

In passing, let us mention that arguments about de- 
velopment of the most optimum effective description of 
hadronic/nuclear physics in terms of mesonic degrees of 
freedom, known also as Quantum Hadro Dynamics, are 
still going on. Using some field variables can be better 
than others: in this respect let us mention the paper 
[8] which emphasized the instead of the traditional a of 
the linear sigma-model, a chiral partner of the pion, one 
better use the radial field \J a 2 + n 2 which has normal 
derivative coupling to pions. There is extensive litera- 
ture on loop corrections and related observables, such as 
resonance mass and shape modification in nuclear and 
excited hadronic matter, see e.g. [9] as a recent example. 

For the purposes of this work it would be sufficient to 
use simple and widely used Walecka model, which keeps 
only the sigma and the omega exchanges in the effective 
nuclear forces 
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B. Strange valence and strange sea of the a meson 

Scalar and pscudoscalar mesons are different from 
more familiar vector and axial ones in terms of their fla- 
vor composition. In the latter case the so called Zweig 
rule applies, forbidding flavor mixing: so for example the 
u> meson we will discuss in his work has a truly negligible 
mixing with a strange counterpart, <j> meson. Scalar and 
pseudoscalar mesons are on the contrary nearly ideally 
SU(3) octets and singlets, so different flavors are very 
strongly mixed together. The pscudoscalar channel has 
been studied extensively, and we know that in this case 
rf is much heavier than (twice more strange) n meson. 
This is the famous Weinberg U(l) problem, resolved by 
existence of the instanton-induced repulsive interaction 
pushing the singlet state upward. 

The same instanton-induced interaction has the oppo- 
site sign in the SU{3) singlet scalar channel, pushing mo- 
downward, see [6] for details. The magnitude of flavor 
mixing matrix element in the scalar channel has been 
also evaluated on the lattice, and the results agree with 
the instanton-based predictions in sign and magnitude. 

Apart of theoretical motivation, there is a simple phe- 
nomenological fact that no purely strange ss counterpart 
to sigma resonance /o (600) is seen. There are strong evi- 
dence that the pair of states /o(980), a (980) are the KK 
molecule, so the next /o are at 1300 and 1500 MeV. 

For those reasons, we think that a description of a as 
a SU (3) singlet state 



a = -^=(uu + dd + ss) 
v3 



(4) 



strongly split from the SU(3) octet one is a reasonable 
approximation. This means that the valence contribution 
to the derivative is 



dm. 



\val 



< a\ss\a >= 2/3 



(5) 
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since with probability 1/3 there is a strange pair. A 
mixing between different scalar mesons fo (<r = /o(600), 
/o(980), /o(1370)) would further change the valence con- 
tributions downward. However, based on large gap be- 
tween sigma and other states, and also based on bet- 
ter studied 7777' mixing, one might think that the relative 
change due to mixing with next fo states is not signifi- 
cant. 

Let us now consider the contribution of the so called 
strange sea, virtual ss pairs, which are always present 
even in a completely non-strange hadrons like a nucleon. 
The relatively well studied case is the nucleon mass sensi- 
tivity to m s , already discussed in [3]. Let briefly remind 
that KN scattering data imply that 



dm 



N 



<9m. 



=< N\ss\N >w 1.5 



(6) 



dlJln 



dm s 



1 



(7) 



As we will see later, the exact value of the common sea 
contribution is not actually important, since its contri- 
butions to a and to mesons tend to cancel each other 
nearly exactly when we calculate variation of N — N in- 
teraction. What matters is the difference between their 
strange seas, to be discussed in the next section. 

It is convenient to present the effect of the possible 
quark mass variation on the a mass in the following form 



to., 
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m u + m d Sm q 

TO., TO„ 



) = 0.4 



(5m, 



m. 



0.04^ 

TOo 



Sttin „,„<5to s nn ,„Sm 
= 0.19 + 0.045 q - 

rriN m s m q 



(9) 



C. KK,nrj mixing with a 

As the second approximation, we will discuss loop ef- 
fects, or mixing with 2-meson states. Those are also 
completely different for a and u mesons. As we men- 
tioned earlier, the latter practically does not mix with 



KK states, while a does mix with them strongly. An 
admixture of virtual KK, 7777 pairs can be viewed as an 
additional contribution to the strange sea, on top of the 
strange content of the non-strange constituent quarks (7). 

The g mixing with continuum of pseudoscalars is de- 
scribed by the standard mass operator given by the usual 
loop diagram 



£(Q 2 



\ 2 (k,Q) 



2tt 4 [(k + Q/2) 2 - m 2 + ie] [(k - Q/2) 2 - m 2 + ie] 



and thus about 1/5 of the nucleon mass comes from the 
strange sea (m s w 120 MeV). 

Similar matrix elements for a, uj mesons is not possible 
to obtain experimentally, although it can be done on the 
lattice. To estimate it we will adopt a simple constituent 
quark picture, assuming additivity of the strange sea. If 
so, the derivatives analogous to (6) for all mesons should 
be 2/3 of it, or 



where A is the KKa, nr/a or wira couplings. Its real part 
describes the shift of m a due to repulsion from these 
states: the imaginary part (for pions only) gives the 
width. The sign of the shift due to KK is obviously 
negative, since m a < 2m^. The effect of tttt has con- 
tributions of both signs. For constant coupling the shift 
is logarithmically divergent, in reality it has to be reg- 
ulated by formfactors in the vertexes. The total shift 
is negative and large, of the order of very large width 
r CT <~ 300 MeV. Note that this large negative shift is 
partly the reason why the a mass is so small. However, in 
this paper we focus on the dependence on quark masses. 
Assuming that the main dependence comes from masses 
of Goldstone bosons, m v ,m K , we differentiate the mass 
operator over these masses and obtain the convergent re- 
sult. Thus one can ignore formfactors and extract the 
effective coupling constant out of the integral. 

For the derivative at Q — m a = 0.5 GeV we get the 
following numerical value for the shifts 



as 



0.0229GeV^ 2 A 2 



oKK 



«9E 

9m 2 



0.019Gey~ 2 A 



2 \2 



(11) 



(12) 



The couplings are not experimentally known, so we rely 
on the SU(3) symmetry and relate them to X a7T7T , which 
is in turn related to sigma meson width 



(8) 



3 A 2 



4to 2 



where we have used m s = 120 MeV. We see that the 
relative change of the strange quark mass produces much 
larger effect than the relative change of the light quark 
mass. This is similar to the case of the nucleon mass 
variation 



2 167rm 



(13) 



Taking r ff =250 MeV we obtain A 2 7r7r = 5 Gev 2 . The fac- 
tor 3/2 account for tt + tt~, 7r°7r° modes. However in the 
mass shift there are contribution of K + K~ ,K K 7 r/ri 
channels which we would count as 5/2. Substituting 
numbers and using standard Gell-Mann-Oaks-Renner ex- 
pressions for m v , rriK (m 2 oc ttIsAqcd), we obtain the ad- 
ditional sensitivity of the sigma mass shift arising from 
the mixing effects 



8m a Sm s 



TOo 



777,, 



TO 



K 



J 2m 2 dm 2 K 



+ 0.5 



rn. 



<9£ 



2m 2 <9m 2 



<5m. 



m. 



0.14 (14) 



We included K + K~,K°K° modes with the coefficient 1 
and 7777 with 1/2: the latter only contribute about 1/5 of 
the final answer. We ignored even smaller contribution 
of the rj'rj' loop. 
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D. The total sensitivity of m CT as compared to mjv 

Together with the one estimated in the previous sec- 
tion, it leads to total 



^« (0.24 + 0.16 + 0.14)^ 



0.54 



5m s 



m s 



(15) 



where 3 terms are the contributions of the common 
strange sea (7), valence strangeness (5) and the loop mix- 
ing (14), respectively. 

In the same units the sensitivity of the nucleon mass 

is 



5m 



N 



m N 



0.19 



Sm s 



m s 



(16) 



We conclude that the sigma mass is about 3 times more 
sensitive to the variation of the strange quark mass than 
the nucleon mass. 

We also need the sensitivity of u> meson to the strange 
mass variation. This meson does not have valence strange 
quarks and practically does not mix with </>, K and rj 
mesons. Therefore, only the strange see contributes: 



0.15- 



. 5m, 



(17) 



III. THE MODIFICATION OF THE DEUTERON 
BINDING 

A. Preliminary analytic estimates 

Simple analytic estimate for sensitivity of the deuteron 
binding to sigma and omega mass modification is ob- 
tained by the differentiation of the potential over the 
mass and averaging the resulting expression <~ exp(-mr) 
over the radial wave function. 

The simplest short range approximation leads to the 
free motion wave function 



ip(r) 



/2k 



-exp(-Kr) 



(18) 



This wave function tends to infinity at small distances. 
However, the real deuteron wave function should be small 
there because of the repulsive core in the potential V(r). 
Therefore, we introduced a small cut-off radius b in the 
integration over r. We estimated b from the condition 
V(b) = which gives b — 0.45 fm. We get the following 
shift of the deuteron binding 



5Qd _ m a Sm a g 2 2k 
Qd Qd m a 4ir2K + m r7 



A variation of m w gives 



5m u 



2k, 



SQd = 

Qd Qd Air 2k + m u 



-(2K+m (7 )b 



-(2K+m„)fc 



-75^ (19) 



m„ 



80- 



5m u 



(20) 



The sign difference between these two derivatives and 
very large derivative value is due to fine tuning between 
omega and sigma terms, which are separately much larger 
than the sum. 

The next step one can do analytically is to add a sim- 
plest square potential to the hard core. The energy of a 
shallow level in such potential is equal to [11]: 



Ed = -Qd = 



Un 



7T 2 (U - Up) 2 

16 U 



8ma 2 ' 



(21) 



(22) 



Here U and a are the depth and width of the potential 
well (a = c — 6, where c and b are outer and inner radii), 
m = m^/2 is the reduced mass. Selecting the width and 
depth of the well to be a = 1.6 fm, U = 40.2 MeV, U = 
52.6 MeV, we got 



SQd 
Qd 



-81.6e 



5m a 



m„ 



^ w 87.4e- m " b 

Qd 



6m u 



m u 



(23) 



(24) 



By changing the core radius from b — to b = 0.4 fm 
one can vary the answer by about factor 3. Simple expo- 
nential dependence on the core radius b appears because 
of translational invariance of Id Scrcdingcr equation for 
rip(r). 

Another effect one should consider is the modification 
of the nucleon mass: its contribution to modification of 
the deuteron binding is 



SQd 

which leads to 
5Q 
Q 



SM 



N 



M N 



< d\ 



jr 



2M N 



\d> 



U + U 5m N 7 Sm N 



U — Uq ttin 



m N 



(25) 



(26) 



Although the sensitivity to the nucleon mass is much 
weaker than that for mesons, it is still quite strong: we 
attribute it to the fact that the small deuteron binding 
energy is in turn the delicate balance between larger ki- 
netic and the potential energies. 



B. Using the Walecka potential 

The estimates of the preceding section are given for 
orientation only, and in fact one of course have to solve 
numerically the radial Schreodinger equation and obtain 
the correct wave function. The one can either average 
the potential derivative over it or simply vary all masses 
involved explicitly. We did the latter and determined the 
sensitivity of Qd to the sigma, omega and nucleon mass. 
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Strictly speaking, at this point it is no longer possible 
to limit ourselves to Walecka model, with the coupling 
constants (3), since it does not describe correctly the 
deuteron binding. In fact, by ignoring all spin-dependent 
forces one cannot even separate the spin-singlet and the 
spin-triplet states. The tensor forces, attributed to pion 
and rho exchanges, are needed for this task. Instead of 
doing so, we have chosen to modify a bit the strength of 
the omega term, reducing g 2 by a factor 0.953 as com- 
pared to (3) and obtaining the correct deuteron binding. 

Our results are* 



SQd ^ 
Qd ~ 

SQd 
Qd 

SQd 
Qd 



-48- 



5m„ 



-26 



5m s 



50 



m„ 



Sm u 



rn. 



7.5 



Sm s 



6 



77V 



5m n 
m N 



1.1 



771, 



Sm s 



m s 



(27) 



(28) 



(29) 



One can see that the first two derivatives are more sen- 
sitive to exact shape of the wave function: they agree 
qualitatively but not quantitatively with the analytic es- 
timates above. We will not show here such details as ex- 
act and approximate wave functions, but just comment 
that the difference between those explain the difference 
in the integrals. 

Summing all the contributions we find 



SQd 
Qd 



-17- 



5m s 



(30) 



Using limits on Big Bang Nucleosynthesis from [3] 
,SQ d 



< 0.1 



(31) 



zero [12]. Today the lowest resonance energy is only 
E = 0.0973 ± 0.0002 eV is large compared to its width, 
so the neutron capture cross section a <~ l/^o- The 
data constrain the ratio of this cross section to a non- 
resonance one (which was used to measure number of 
neutrons emitted by the reactor). It therefore implies ^ 
that these data constrain the variation of the following 
ratio 5{Eq/E\) where E\ <~ lMeV is a typical single- 
particle energy scale, which may be viewed as the energy 
of some 1-body "doorway" state. 

A generic expression for the level energy in terms of 
fundamental parameters of QCD can be written as fol- 
lows 

SE — A * SKqcd + B q * 5m q + B s * Sm s + CSa * A QC d (33) 

where A, B, C are some coefficients. The first term is 
the basic QCD term, while others are corrections due to 
modification of the quark masses and the electromagnetic 
a. 

In this section we provide new estimate of the B s . 
More specifically, we estimate the variation of the res- 
onance energy resulting from a modification of the sigma 
mass. The energy of the resonance E = E excita tion ~ S n 
consists of excitation energy of a compound nucleus, mi- 
nus the neutron separation energy S n . This, in turn, is 
a depth of the potential well V minus the neutron Fermi 
energy cf, S n = V — 6f- The latter scales like 1/R 2 if 
the radius of the well is changed. The kinetic part of the 
excitation energy E excit ation scales in the same way. So, 



Eq — E exc itation S n E, 



excitation 



+ e F -V = K 



MR 



V (34) 



where K is a numerical constant which can be found from 
the present time condition Eq ~ 0. The shift of the 
resonance then is 



one gets the final limit on the m s variation to be 
5(m s /k QCD ) 



(m s /AQCD) 



< 0.006 



(32) 



ft 2 (SM 2SR\ CTA / SM 2SR SV 

SE = -K-— — + —— - 5V = -V — + — - + 



MR \ M R 



V M 



R 



Using eq. (2) we can find the depth of the potential well 



V 



IV. OKLO 

In this section we extract limits on Sm s following from 
data on natural nuclear reactor in Oklo active about 2 
bn years ago. The most sensitive phenomenon (used 
previously for limits on the variations of the electro- 
magnetic a) is disappearance of certain isotopes (espe- 
cially Sm 149 ) possessing a neutron resonance close to 



*In the derivative over the omega mass we have divided back 
by the factor 0.953, restoring the original strength of the vec- 
tor term. 



V = 



in- 



A 

m 2 , 



(36) 



Here ro~ 1.2 fm is an inter-nucleon distance. Numerical 
estimates shows that the contribution of the variation of 
r (and the variation of R = A 1 ^ 3 ^ in eq. (35)) is not as 



*Of course, under assumption that the same resonance was 
the lowest one at the time of Oklo reactor. 

*Note that the suppression of N-N wave function at small 
separation due to the repulsive core reduces the depth of the 
effective potential V. However, this effect is not so important 
in the ratio 
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important as the direct contribution of the m a variation 
in the equation above. This gives us 



SV „ „dm a „ „6m u 
— = -8.6 — - + 6.6 — - 
V m a m w 



= -3.5 



Sm s 



m s 



SE Q = 1.7 x 10 8 eV x 



5m. 



(37) 



(38) 



We used V= 50 McV in cq. (35). Comparison of this 
result with the the observational limits claimed in [12] 
\6Eq\ < 0.02 eV gives a very strong limit 



Mm s / A QC p) 

(m s AQCD) 



< 1.2 x 10 



-10 



(39) 



at time w 1.8 billion years ago. 

Note that the authors of the last work in [12] found 
also the non-zero solution 5E a = -0.097 ±0.008 eV. This 
solution corresponds to the same resonance moved below 
thermal neutron energy. In this case 



S(m s /A Q cp) 
(m s A QC D) 



= -(0.56 ±0.05) x 10" 



The error here does not include the theoretical uncer- 
tainty. 

The production of nuclei with A > 5 during BBN is 
strongly suppressed because of the absence of stable nu- 
clei with A = 5. 5 He is unstable nucleus which is seen as a 
resonance in n-a elastic scattering. The ground state lies 
at 0.89 MeV above neutron threshold. The variation of 
the strange quark mass may influence the position of the 
resonance making, for example, 5 He stable. Stable 5 He 
at the time of BBN would change strongly the primor- 
dial abundances of light elements. The estimate similar 
to that we made for Sm nucleus gives us a limit 



S(m s /A Q cp) 
{m 8 A QC D) 



> -0.006 



(41) 



This limit corresponds to SE > - 0.89 MeV at the time 
of BBN. 

We obtained limits on variation of m s /AQCD during 
the interval between the Big Bang and present time and 
on shorter time scale from Oklo natural nuclear reactor 
which was active 1.8 billion years ago. It is also possi- 
ble to obtain limits on the intermediate time scale. One 
possibility is related to position of the resonance in 12 C 
during production of this clement in stars. This famous 
resonance at E=380 KcV is needed to produce enough 
carbon and create life. According to Ref. [13] the posi- 
tion of this resonance can not shift by more than 60 KeV 
( one can also find in Ref. [13] the limits on the strong in- 
teractions and other relevant references). We have made 
a very rough estimate of the limit on the strange quark 
mass variation which can be obtained via m^-mcchanism: 



8{m s /A QCD ) 

(m 8 AQCD) 



< 0.001 



(42) 



This limit can be improved after an accurate calculation. 
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